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1 . INTRODUCTION 

Over the past several years we have conducted many experiments In the 
area of global satellite data assimilation. Most of these experiments have 
dealt with the assimilation of satellite-derived temperature data either 
continuously or Intermittently Into a global primitive equations model. We 
have found that when properly performed, satellite data ass lml la cl on has a 
positive Impact upon the global model forecasts. Using the experience we have 
gained from our global satellite data assimilation experiments, we have turned 
our attention to the problems of mesoscale satellite data assimilation. Our 
overall research goal is to utilize the assimilation of satellite data in 
order to improve the forecasts made by mesoscale forecast models. 

Specifically, the problem we report on here has dealt with the 
assimilation of high resolution (~30 km) satellite-derived temperature data 
Into a mesoscale model with horizontal resolution of 50-60 km. Unlike global 
assimilation. In which a small portion of the forecast model domain Is subject 
to data Insertion at virtually every time step, mesoscale assimilation Is such 
that virtually all of the forecast model domain is subject to data Insertion 
at one time step. Thus, the mesoscale problem lends Itself natually to 
intermittent data assimilation and the forecast model is re-lnitialized 
whenever a new satellite pass covers Its domain with data. In this research 
we have treated the satellite data assimilation as an initialization problem. 

In the next section we discuss the satellite and conventional data sets 
available for use in this study. The variational temperature analysis scheme 
we have developed to combine the temperature Information contained in these 
two types of observations Is described and demonstrated in Section 3. The 
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CIMMS Mesoscale Model Is the forecast model which we used In this study to 
assess the Impact of the satellite data assimilation. It evolved from the 
global primitive equations model we have used in our global assimilation 
work. A description of the model is given in Section 4. Finally, in Section 
5 we discuss the results of the forecast experiments which culminated our 
research efforts. Three different initialization procedures are compared with 
and without satellite data assimilation. 


/ 


I 


-3- 




2 . DATA ASSESSMENT 
Satellite Data 

The satellite data which were available for this study were brightness 
temperatures from one orbit for the TIROS-N High Resolution Infrared Radiation 
Sounder (HIRS) and Microwave Sounder Unit (MSU) channels as determined by NESS 
at the University of Wisconsin. The TIROS-N orbit was for 15 January 1979 at 
0850 GMT and covered roughly the eastern half of the United States. The 
resolution of the satellite data was on the order of about 30 km. 

In Fig. 1, It can be seen that at 12 GMT, 15 January 1979 a large high 
pressure system was centered on the eastern U.S. Ab a result, moBt of this 
area was relatively cloud-free as indicated by the nephanalysls displayed In 
Fig. 2. The density of the TIROS-N observations Is compared to that for the 
FGGE Level Il-b surface observations In Fig. 3. The satellite observations 
shown in this figure are those from areas that have been determined to be 
relatively cloud-free. In this study we are Interested in utilizing the 
Information about horizontal temperature gradients contained In the TIROS-N 
brightness temperatures. From Fig. 3 we can see that one can expect gradients 
to be more sharply defined by the higher resolution satellite observations 
than the surface observations. The upper air observational network has even 
lower resolution than the surface network. 

We next examined the brightness temperatures for the various HIRS and MSU 
channels. Plot6 of this variable for HIRS Channel 13 and MSU Channel 3 are 
shown in Figs. 4 and 5. The brightness temperature for HIRS Channel 13 can be 
thought of as roughly a mean layer temperature for a layer extending from the 






Fig. 2. 


Observed sky cover and nephanalysis from surface observations 
for 12 GMT, 15 January 1979. 
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Fig. 4 


HIRS Channel 13 brightness temperatures 
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eurface to about 800 mb while that for MSU Channel 3 represents the mean layer 
temperature for a layer extending from about 600 mb to 100 mb. These infer- 
ences are based upon the weighting functions for these channels shown in Fig. 
6. In the cloud-free areas we can see in Fig. 4 that temperature gradient 
structure is evident while in Fig. 5 we see that there is virtually no temper- 
ature gradient whatsoever. This can be understood by examining the weighting 
function for MSU Channel 3 and the FGGE Level Ill-b temperature fields valid 
at 12 GMT, 15 January, at 500 mb, 300 mb, and 150 mb displayed in Figs. 7-9. 
In Fig. 7 we see that there is a definite north-south temperature gradient 
with cold air to the north at 500 mb. This gradient is still apparent at 300 
mb though not as strong (Fig. 8), but at ' *1 mb we see in Fig. 9 that the 
temperature gradient has reversed itself with warm air to the northeast and 
cold air to the southwest. The broad weighting function for MSU Channel 3 
indicates that the resulting brightness temperatures will represent an average 
for a deep layer. The horizontal temperature gradients at the top of this 
layer are almost the reverse of those at the bottom of the layer. Thus, the 
resulting average is a field with virtually no horizontal temperature 
gradient. Therefore, we concluded that the TIR0S-N brightness temperatures 
themselves were of most use for the very lowest layers of the atmosphere where 
the weighting functions for the appropriate channels peak more sharply. 

FGGE Level Il-b and Ill-b Data 

In thi 8 study we had available both FGGE Level Il-b and Ill-b data for 
the time period of interest. The FGGE Level Il-b data set is the Global 
Experiment Research Data Set, which is distinguished from the operational Il-a 
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Fig. 6 


Weighting functions, for U.S. standard atmosphere at nadir viewing 
for seven HIRS2-MSU channels, (from Sussklnd, et al., 1982) 
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FGGE Level Ill-b 300 mb temperature field, 12 GMT, 15 January 1979. 
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datn set by a delayed cutoff to acquire a complete global data set. This data 
aet contains the majority of all routine weather observations globally 
observed and transmitted over the Global Telecommunications System In 6-hour 
cycles. These data are augmented on a massive scale with observations from 
many special systems: satellites, aircraft, buoys, ships, and balloons during 
the FGGE Special Observing Periods. The FGGE Level Ill-b data aet consists of 
analyses produced by the European Centre for Medium Range Weather Forecasts 
operational system using the Level Il-b rather than Il-a data sets as Input. 

Using the FGGE Level Il-b and Ill-b data sets from 12 GMT, 15 January 
1979, we analyzed temperature at the surface as well as at the standard pres- 
sure levels (1000 mb, ftiO mb, 700 mb, etc.). Heights and temperatures at the 
standard pressure levels were obtained from the FGGE Level IXI-b data, and 
using the terrain hel/hts Illustrated In Fig, 10, a surface temperature field 
was produced by Interpolation to a 0.5° grid for the domain Illustrated. This 
surface temperature field was then used as the first-guess field for a 
Cressman analysis of the Level Il-b surface temperature observations. The 
resulting temperature analysis Is shown In Fig. 11. The radii of Influence 
for the Cressman analysis were chosen to reflect the resolution of the surface 
observations as shown In Fig. 3. We can see from this temperature analysis 
that the strongest temperature gradients are found along the coastal areas of 
the Gulf of Mexico and the southeast Atlantic seaboard. The definition of 
these surface temperature gradients is limited by the resolution of the sur- 
face observations. Using the FGGE Level Ill-b temperature fields a6 first- 
guess fields, Cressman analyses of the Level Il-b upper air observations were 
performed for the standard pressure levels. The radii of influence for these 



Fig. 10 Terrain height used for surface analyses. 




analyses were chosen to reflect the resolution of the upper air data. Figs. 
12 and 13 display the resulting analyses for 700 mb and 200 mb. 

In view of our examination of the satellite data we concluded that It 
would only be feasible to produce temperature analyses enhanced by TIROS-N 
brightness temperatures in the very lowest layers of the atmosphere. Using 
the analyzed surface temperature field shown in Fig. 11 as a first-guess 
field, the brightness temperatures from HIRS Channel 8 (the 11 pm longwave 
window) were analyzed. Only the brightness temperatures from areas determined 
to be relatively cloud-free were used (see Fig. 3). The resulting Cressman 
analysis, whose radii of influence reflected the high resolution of the satel- 
lite observations, is shown in Fig. 14. The most striking difference between 
this satellite-derived surface temperature analysis and the conventional 
analysis shown in Fig. 11 is the enhanced temperature gradients in the coastal 
areas. These gradients are much more sharply defined in both their horizontal 
structure as well as their intensity. This is graphically illustrated in Fig. 
15 where the difference field between these two analyzed fields is shown. 

In order to better understand the differences between these two surface 
temperature analyses, each was compared with the actual surface temperature 
observations. The bias and rmse were computed for the differences between the 
surface observations and the values of each field interpolated to the loca- 
tions of the observations. For the conventional analysis the bias and rmse 
were -0.01K and 1.04K, respectively. The respective values for the satellite- 
derived analysis were -1.58K and 3.76K. Since the first field is merely a 
Cressman analysis of the surface observations themselves, we expect that the 
values for the bias and rmse will be quite small, and indeed they are. The 















Fig. 13 Cressman analysis of FOGE Level II-b upper air 200 mb temperature 
observations, 12 GMT, 15 January 1979. 
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satelll te-deri ved field displays a large bias with respect to the surface 
observations along with a larger rmse. However, the high resolution tempera- 
ture gradient Information contained in this field is potentially valuable. In 
the next section we describe a variational analysis scheme designed to combine 
the useful information contained In the conventional and satellite-derived 


surface temperature fields. 


3. VARIATIONAL TEMPERATURE ANALYSIS 


The philosophy behind the variational temperature analysis scheme la to 
combine horizontal temperature gradient Information derived at high resolution 
from satellite with absolute temperature Information derived from conventional 
observations and/or a model forecast. The analysis Is obtained by minimizing 
the functional J in spherical coordinates: 

0 A 

J - f 2 / 2 f A(T-T) 2 + B(V.T-V.T) 2 + C(7 n T-V fl T) 2 

0 X, A A 06 

1 1 ( 1 ) 

+ D(V 2 T-7 2 T) 2 + E(V 2 T-V?T) 2 } dAde 
A a o y 

where T Is the analyzed temperature, T is the absolute temperature, and T is a 
satellite-derived temperature. The gradients and Laplacians are separated 
into their longitudinal (X) and latitudinal (0) components. A, B, C, D, and E 
are data weights whose values are chosen to be roughly equal to the inverse of 
the expected variance of their associated differences in (1). The units of 
these weights are such that each term in (1) Is dimensionless. In order to 
satisfy the boundary conditions when the variation of J is taken, the weights 
are set to zero as appropriate. This functional can be evaluated on any 
surface and is readily solved by relaxation. 

In order to test the analysis scheme we use it to combine the surface 
temperature field derived using conventional observations displayed in Fig. 11 
with the field shown in Fig. 14 derived using high-resolution satellite obser- 
vations. The conventionally analyzed temperatures are denoted by T and the 

A 

satellite-derived temperatures by T in functional J. The variational 


temperature analysis Is to be performed upon a 41 x 51 grid with 0.5° 
resolution extending from 25 N to 45 N and from 95 W to 70 W, 

By varying the relative values of the weights In (1), we can control the 

* 

Influence of the T- and T-flelds upon the analysis resulting from minimizing 
the functional J. Three cases will be shown here. In each case the value of 
A Is taken to be 0.1. In the first case B - C - 1.E9 and D - E - 1.E18. The 
resulting temperature analysis (Tj) Is shown In Fig. 16 and the difference 
field (Tj-T) In Fig. 17. Table 1 summarizes some statistics computed for the 
three cases (Tj, T 2 , and T^) as well as for the T- and T-flelds. In the first 
case the relative weighting upon the satellite temperature gradient and 
Laplaclan Information Is weak and the analysis field agrees more with the 
surface observations. Cases 2 and 3 are shown In FigB, 18-21 and represent 
medium and strong weighting of the satellite temperature Information, respec- 
tively. In Case 2, B ■ C * 3.E9 and D ■ E - 3.E18 

Table 1. RMSE's between various temperature fields and observations of 
temperature and horizontal derivatives of temperature. 


Field 

Surface Obe. (K) 

VT(K» _1 ) 

V T(K* £ ) 

a/ 

T 

1 .04 

2.7E-5 

t-t . IE-10 

T 

3.76 

0 

0 

T 1 

1.70 

1.2E-5 

1.3E-10 

T 2 

2.23 

7.1E-6 

6.0E-11 

T 3 

2.70 

3.6E-6 

2.5E-11 


while In Case 3, B - C - 1.E10 and D « E > 1.E19. From the RMSE's In Table 1 
we can see that as the relative strength of the weighting upon the satellite 





Fig. 16 


Variational temperature analysis (Tj ) for 12 GMT, 15 January 1979 
(A - 0.1, B - C - 1.E9, D - E - 1.E18). 





Fig. 17 


Difference field (Tj - T) betw 
15 January 1979. 


analyzed temperatures for 12 GMT 













Fig. 19 Difference field (T 2 ~ T) between analyzed temperatures for 12 GMT 
15 January 1979. 




Fig. 20 


Variational temperature ana 
(A « 0.1, B - C « 1.E10, D 






temperature information increaaee, the resulting fields are such that their 
gradients and Laplacians more closely match those for T-field while at the 
same time they less closely match the surface observations of temperature. As 
the difference field in Fig, 15 would suggest, it can be seen in Figs. 17, 19 , 

and 21 that the major impact of the satellite information is upon the coastal 
areas of the southeastern U.S. 
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4 . ME SOS GALE MODEL 

The forecast model used In thia study is a dry version of the CIMMS 
Mesoscale Model (CMM). It is a limited-area primitive equation model (with 
terrain) which evolved from the global model described by Sasaki and Goerss 
(1982). A unique feature of this model is that it employs variational con- 
straints each time step in order to control the total mass and energy within 
Its domain (Sasaki, 1976; 1977). In this section we will briefly outline the 
features of the CMM. 

The model uses spherical coordinates with a as the vertical coordinate. 
It employs a staggered latitude-longitude grid using the Arakawa (C) lattice 
(Mesinger and Arakawa, 1976) as illustrated in Fig. 22 for 0.5° grid spacing. 
Values of tt , T, ft, and $ are carried at the grid poinls labeled A. The vert- 
ical structure of the model is also shown in Fig. 22. The number of levels 
where v and T are carried 1b denoted by KM. The value of a is zero at the top 
of the model atmosphere (p t ) and is one at the surface (p g ). The vertical 
velocity 6 is zero at both a - 0 and o ■ 1. The model consists of KM layers 
each 1 /KM thick. 

The primitive equations are written in flux form using (X, 0, o) 
coordinates : 


3n + 1 

3t "a cos 0 


{gy(iTu) + 
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V cos 0)} + 
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3t 


Uu) + a'coiTe { h (lru2) + fe ( " uv cos + lo (l,u5) 


n .3* , 3n. , u tan 0 V _ 

- — rrr + era—) - (f + — ) ttv - ttF, 

£t cos 0 oA 9 A a \ 
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Fig. 22 Horizontal and vertical structure of the CIMMS Mesoscale Model. 
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TF (vV) + a ~cos~ 6 { h ( " UV) + I? ( " v2 C ° S 6)) + h (vvh) 
, tt ,3* . 3^\ , /£ . u tan 6> 

I ( 99 + CTa Te ) + (f + — i — > nu - " F 0 


ft (ttT) + a 'cos e { bt vuT) + fe <tfvT cob 6)} + ( Tooo )K fl ( ™ 6) 


+ ■ tt2 - oq r- 

c '■a cos 0 3A T 90 
P 


1 9 u 3 , 3o i tt Q 

IT + 30< V COfl + 9?] “c 


where n ■ P 8 - P t and p - tid + p fc . The value of p t Is 150 mb. Except for 
friction in the lowest o-layer, no parameterization of irreversible physical 
processes was included in this study. The model integrations are carried out 
using centered -time and centered-space finite differencing in order to fore- 
cast tt , u, v, and T. After the forecast for it at time n+1 has been made, 
diagnostic equations are solved for the values of a and at time n and the 
rejuits are used in forecasting u, v, and T at time n+1. 

The CMM employs only very simple and computationally compact finite 
differencing while variational adjustments are made each time step in order to 
control the total mass and energy within the model domain. The values for 
these quantities are diagnosed either from the larger scale model (if any) the 
CMM is nested within or from the analyzed fields which provide the CMM with 
its lateral boundary conditions. 

For a given time step, in order to control the mass within the CMM domain 
we minimize the functional; 

(ff-Tr)^A/, + A m (J \ ttAA - M), 
i j ij 
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where M Is Che total mass within the model domain at that time step as deter- 
mined from the larger scale model or the analyzed fields, AA Is the Increment 
of area associated with each model grid point, tt Is the current model solu- 
tion, and ti Is the adjusted value. The resulting correction to the 71 -fleld 
Is : 


71 ■ IT - X m /2 

where X - 2(£ ) ttAA - M)/ ][ J AA. 

1 j i j 

This correction Is applied at each time step. 

To control the total energy within the model domain each time step, the 
following functional is minimized: 


J ■ [ | I MCu - ^) 2 AA u + (v - v) 2 AA v ] + 0(T - T) 2 AA^,} 


+ X-. [J T y (c u^TAA^ + ti u 2 AA /2 + tt v 2 AA /2) 
E 'J j J ' p I T u u v v 


- TEi , 


where TE Is the total energy within the model domain at that time step as 
determined from the larger scale model or the analyzed fields; AA Is the 
Increment of area associated with each model grid point; u, v, and T are the 
current model solutions while ti Is the result of the previous variational 
adjustment; and u, v, and T are the adjusted values. The weights a and 6 are 
chosen so that their ratio Is T/c p ■ 0.272. The resulting corrections for u, 


v, and T are: 






4- v/ c. 
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u - u(l - X £ Ti u /2a) , 

v - v(l - X E n v /2a) , 

and T » T ~ A„c n„/2?> . 

E p i 


An approximate solution 


for A Is used: 
E 


I l l <Vt ¥ &A T + % A V 2 / 2 + % AA v'' 2/2) ‘ TE 

_ 11k 

E 111 <y 2 AA t /28 + AA u u 2 /2o + ti 2 AA y v 2 /2a) 

1 j k 


Finally, In order to compensate for the adjustment to u and Its effect upon 
the energy calculation, X Is adjusted by the ratio £ j ttAA/M before the 

E i J 

adjustments for u, v, and T are computed. 

Normally a "double-nesting" procedure Is employed with the CMM nested 
within a coarser resolution version of itself which in turn can either be 
nested within the global model or have its boundaries specified in Borne other 
fashion. For example, a 1° resolution version of the CMM might be used to 
make a 12-hour forecast beginning at 12 GMT with its boundary tendencies 
determined from analyzed fields at 12 GMT, 18 GMT, and 00 GMT or from the 
forecast fields obtained from a 2° resolution version of the global model. 
Then forecasts could be made using a 0.5° resolution version of the CMM whose 
domain would be in the center of the 1° version's domain and whose boundary 
tendencies are determined from the forecast field produced by the 1° model. 


The lateral boundary conditions employed In the CMM are those developed 
by Perkey and Kreltzberg (1976). They allow changes outside the limited 
domain to Influence results there without contaminating the forecast with 
spurious boundary reflected energy. Their boundary condition transforms long- 
and medium-length Interior advectlve and gravity waves Into short waves that 
can be removed by a low-pass filter. For any prognostic variable X and time 
n+1 , 


3X 


- W> + w 


3X 

I.n « + <1 - "(») jrl I>n « 


where 1 specifies the large-scale tendency of X and m specifies the model 
tendency. The weights are given by; 


W(I) - 


0.0 I * boundary grid points 
0.4 I - (boundary-1) grid points 

0.7 1 - (boundary-2) grid points 

0.9 I ■ (boundary-3) grid points 

1.0 I “ all other points 


The low-pass filtering is accompllsed by using a smoother-desmoother upon 
the grid points within 7 grid intervals of the boundaries every Nth and 2Nth 
time step and over the entire domain every 3Nth time step. The choice of N is 
dependent upon the model resolution and the length of its time step. In one- 
dimension the smoother-desmoother used 1 b: 


xj - (1 “ a)X 1 + | (X 1 _ 1 + X 1+J ) 

where a ■ 0.50 on the first pass and a ■ -0.52 on the second pasB, 

We now present the results for a 12-hour forecast using the dry version 
of the CMM with 0.5° grid spacing and twelve 0 -layers between the surface (p„) 
and 150 mb (p t ). The 0.5° resolution model was nested within a 1.0° resolu- 
tion dry vernion of the CMM with the same vertical structure. The horizontal 
domains of the models are shown in Figs. 23 and 24. The 1.0° resolution model 
was initialized at 12 GMT on 15 January 1979 using fields interpolated from 
the FGGE Level Ill-b data set valid et that time. The FGGE Level Ill-b 
analyses were produced using the European Centre for Medium Range Weather 
Forecasts (ECMWF) operational system and were available globally at 1,675° 
resolution and at the 15 standard levels between 1000 mb and 10 mb, 
Inclusive, The boundary tendencies for the 1.0° resolution CMM forecasts were 
determined by first constructing fields at 18 GMT, 15 January 1979 and at 00 
GMT, 16 January 1979 interpolated from the FGGE Lev Ill-b data sets valid at 
those times. The tendencies were then found by simply interpolating between 
the 12 GMT and 18 GMT fields and between the 18 GM1 and 00 GMT fields. In 
Fig. 23 we see the initial 1000 mb and 500 mb height fields for the 1.0° 
resolution model at 12 GMT, 15 January 1979 along with the forecasts valid at 
18 GMT, 15 January 1979 and 00 GMT, 16 January 1979, respectively. We can see 
in this figure that the Perkey-Krei tzberg boundary scheme appears to be quite 


effective 
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^g. 23 Initial (top), 6—hour forecast (middle), and 12—hour forecast 

(bottom) 1000 mb and 500 mb height fields from 12 GMT, 15 January 

1979 run of the 1.0° resolution, dry version of the CIMMS Mesoscale 
Mod el . 
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Initial (top), 6-hour forecast (middle), and 12-hour forecast 

\ h °l t0 *i JJ°!' mh *f 0 500 “ b bel « ht fl elds from 12 GMT, 15 January 
1979 run of the 0.5 resolution, dry version of the CIMMS Mesoscal 









Fig. 25 


Interpolated 1000 mb and 500 mb height fields from the FGGF Level 
Ill-b data sets valid at 12 GMT, 15 January 1979 (top), 18 GMT, 15 
January 1979 (middle), and 00 GMT, 16 January 1979 (bottom). 
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The 0.5° resolution model was Initialized at 12 GMT, 15 January 1979 with 
fields Interpolated from the 1,0° resolution model. The Initial 1000 mb and 
500 mb height fields are shown In Fig. 24 along with the forecast fields valid 
at 18 GMT and 00 GMT, respectively. The boundary conditions for the 0.5° 
resolution model were obtained by linearly Interpolating between forecast 
fields produced by the. 1.0° resolution model at 2-hour Intervals between 12 
GMT and 00 GMT. Again we see that the Perkey-Krel tzberg boundary scheme works 
well. 

Finally, In Fig. 25 verification fields interpolated from the FGGE Level 
Ill-b data sets valid at 12 GMT and 18 GMT, 15 January 1979 and at 00 GMT, 16 
January 1979 are displayed. Comparing Figs. 24 and 25 we can see that the CMM 
forecasts of 1000 mb and 500 mb heights are quite reasonable. 
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5- ASSIMILATION EXPERIMENT 

Ultimately, one would like to use satellite temperature data to improve 
the forecasts made by mesoscale models. In the previous sections we have 
examined the conventional and satellite data available, a variational temper- 
ature analysis scheme designed to combine these data, and the CMM which will 
be used to test the impact of data assimilation. Since the TIROS-N data is 
such that one can normally expect virtually all of a mesoscale model domain to 
be covered by one pass of the satellite, the data assimilation problem is 
really an initialization problem. Whenever a satellite pass occurs the entire 
model domain would be re-initialized at that time. Thus, we have an intermit- 
tent assimilation proLiera rather than the continuious assimilation problem 
that one normally encounters when dealing with satellite temperature data 
assimilation Into a global model where a small portion of the model domain is 
modified each time-step as new satellite data become available. In this 
section v*e shall describe three different initialization procedures for the 
CMM and test their effectiveness at assimilating the satellite temperature 
information into the CMM. Using the data from 15 January 1979 the impact of 
satellite temperature assimilation will be evaluated by comparing CMM fore- 
casts made with and without satellite data for all three different initiali- 
zation schemes. 

Initialization Procedures 

Each of the model Initialization procedures used in this study consists 
of some form of objective analysis performed upon the standard pressure sur- 
faces at the horizontal grid-points used by the model. The fields used to 
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initlalize Che model are then produced by vertical Interpolation from these 
objectively analyzed fields to the model's o-coordinate system. 

The first Initialization procedure we shall consider Is the three- 
dimensional multivariate variational objective analysis procedure described by 
Sasaki and Goerss (1982). Objectively analyzed fields of geopotential height 
and winds are produced by minimizing the functional 

J " /// [ A( (f>— < f>) 2 + B(~ - T) 2 + C u (u-u) 2 + C y (v-v) 2 

+ d(?$-V$) 2 + e(V 2 <j>-V 2 ?) 2 + f(£-?) 2 + g(D-D) 2 

+ h( a~cos~6 n ~ fv)2 + 1( a fe + fu) ^ dXd6d,r * 
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> V\ *” la ' 


30 
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The resulting values of geopotential height, u-wind, and v-wind are denoted 
by $ , u, and v, respectively, while their corresponding first-guess field 
values are denoted by $ , u, arid v. Three-dimensional data weight matrices A, 
B, C u , and C v contain relatively high values at locations of observed 
geopotential height, temperature, u-wind, and v-wind ($, T, u, and v) and 
negligible values elsewhere. The weights d, e, f, g, h, and 1 are constant 
over most of the domain but are set to zero near the boundaries in order to 
satisfy certain boundary conditions. Ej ch weight is assigned a value equal to 
the inverse of the expected variance o): the difference it is associated with 
in the functional J. 

The solution of this variational problem over the domain results in 
fields of geopotential height, u-wind, and v-wind: 

1) which agree with the surface and upper-air observations; 

2) for which the horizontal gradients and Laplacians of geopoten- 
tial height and the vorticity and divergence of the wind field 
match those for the first-guess fields; 

3) for which the vertical derivatives of the geopotential height 
field are consistent with the observed temperatures; and 

A) which are in approximate geostrophic balance. 

In this study the first-guess fields were interpolated from the FGGE Level 
Ill-b data set valid at 12 GMT, 15 January 1979. Surface and upper-air obser- 
vations from that time were taken from the appropriate FGGE Level 11-b data 
set. Satellite data assimilation was accomplished by replacing the first- 
guess geopot ential height field (<J>) with one derived using the temperature 
fields produced from the variational temperature analysis scheme described in 
Section 3 and by increasing the relative weighting upon the terms in 
functional J containing gradients and Laplacians of $ . 
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The second initialization procedure can be described as a simple two- 
dimensional univariate objective analysis scheme patterned after that 
described by Cressman (1959). For each variable (surface pressure, tempera- 
ture, u-wind, and v-wind) a Cressman analysis is performed using radii of 
influence based upon the average distance between observations. Surface 
analyses are produced for pressure and temperature using smaller radii of 
influence than the analyses of temperature, u-wind, and v-wind performed upon 
the standard pressure surfaces using the upper air observations. These ana- 
lyses are performed one variable at a time and one level at a time with no 
coupling between the variables or the levels. The model is then initialized 
by interpolating these fields to the model's o-coordinate system. 

The third initialization procedure we shall investigate is identical to 
the second except that after the individual fields have been analyzed a static 
balancing procedure is employed. At each level the following variational 
functional is minimized: 

1 * / / [A(<|>-<(J) 2 + B(u-u) 2 + c(v-v) 2 

Vi 

♦ d( rire H- - fv)2 + *<t H + fu > 2 ) 

The values of geopotential height, u-wind, and v-wind derived from the uni- 
variate Cressman analyses are denoted by 4>, u, and v. The resulting values 
($, u, and v) are such that approximate geostrophic balance is achieved at 
each level. Thus, at each level the mass and momentum variables are coupled. 
As we shall see later, this coupling provides the forecast model with fields 


that are In better balance so that the model integrations are Initiated with 
less gravity wave noise. 

Model Forecasts 

A total of six forecast experiments beginning at 12 GMT, 15 January 1979, 
were conducted for all combinations of the three different initialization 
procedures with and without satellite temperature data assimilation. The 
experiments are denoted by two letters with the first letter standing for the 
initialization procedure (V-three-dlmenslonal , multivariate variational; C- 
uni variate Cressman; C-uni variate Cressman with static balancing) and the 
second letter indicating whether satellite data assimilation has been per- 
formed (S-ratelllte data assimilation; N-no assimilation). Before discussing 
the results of the forecapt experiments we shall first take a look at the 
effectiveness of the three Initialization procedures in terms of providing the 
model with balanced fields. 

An effective measure of the gravity wave noise within a numerical fore- 
cast model is the root mean square pressure tendency determined over the model 
domain time step. Ab described by Miyakoda et si. (1978) a major problem of 
data assimilation and of model Initialization Is the balance between the 
model's mass and momentum fields. When a model is initialized such that these 
fields are not in geostrophic balance, a discontinuity is created with the 
current modal solution. The dynamical characteristics of the model fluid 
enable it to remove this discontinuity in such a way as to maintain a state of 
approximate geostrophic balance. This geostrophic adjustment process results 
in the generation and dispersal of inertial gravity waves. Ideally one would 


like to Initialize a numerical model with mass and momentum fields which not 


only match the observed values of temperature, pressure and wind but also are 
balanced to the extent that the model can begin Integration with a minimal 
amount of noise. In Fig. 26 are plotted the root mean square pressure 
tendencies for the three assimilation forecasts (VC,CS, and GS) from zero to 
twelve hours of forecast time. We can see that when the forecast model Is 
Initialized using the univariate Cressman procedure the level of gravity wave 
noise starts considerably higher than for the other two procedures. After 
about 4-5 hours of model integration the noise level has reached that of the 
other methods. While the variational method haa an Initial noise level higher 
than the Cressman procedure with Btatlc balancing, it quickly reaches compar- 
able levels and after about four hours of Integration reaches a level slighcly 
less than the two Cressman techniques. Thus, as we would expect, the methods 
which provide for the coupling of the initial mass and momentum fields produce 
forecasts with reduced levels of gravity wave noise and perhaps the lowest 
noise levels are realized by the three-dimensional variational initialization 
procedure which also provides coupling between the vertical levels of the 
model atmosphere. 

Now we shall begin discussion of the forecast experiments themselves. In 
every case the CMM was Initialized at 12 GMT, 15 January 1979, using the 
surface and upper air observations available from the FGGE Level Il-b data set 
and using as first-guess fields Interpolated fields from the FGGE Level Ill-b 
data set. For all Initialization procedures objectively analyzed fields were 
produced at the following pressure levels 1000, 925, 850, 775, 700, 600, 500, 
400, 350, 300, 250, 225, 200, 175, and 150. The horizontal grid points used 


RMS Pressure Tendency (mb /sec) 







^w#WW«W •« 


-50- 


by each scheme were the same as Is used by the CMM. For the assimilation 
cases the following was done. At the two lowest pressure levels (1000 mb and 
925 mb) the Cressman temperature analyses for experiments CS and GS were 
replaced by the temperature analyses produced by the variational temperature 
analysis procedure. As discussed previously for the variational case (experi- 
ments VS and VN) the first-guess geopotential height field ($) was replaced by 
one derived using the temperature fields at 1000 mb and 925 mb produced by the 
variational temperature analysis scheme. For the non-assimilation experiments 
(VN, CN and GN) all initial fields were produced by utilizing the respective 
initialization procedures with conventional observations alone. 

We shall now examine the initial fields produced by the different proce- 
dures. Since assimilation is only performed in the two lowest levels (1000 mb 
and 925 mb), the only real differences result there. In Figs. 27-32 are 
plotted the 1000 mb height and temperature fields for experiments VS, VN, GS, 
GN, CS, and CN, respectively. In each case we can see the differences between 
the initial tenperature fields for the assimilation and non-assimilation 
experiments. The variational temperature analysis scheme, as we saw in 
Section 3, produces low level temperature fields with enhanced gradients in 
the coastal areas of the southeastern U.S. Comparing the Cressman techniques 
with the variational one we see that both CresBman techniques produce a deeper 
closed low in the 1000 mb height field off the Florida Coast while the vari- 
ational techniques produce a weaker inverted trough. In each case compare the 
differences between the initialized fields and the first-guess fields inter- 
polated from the FGGE Level Ill-b data set shown in Fig. 33. As we can see 
the first-guess fields contain very little detail for either the 1000 mb 
height or temperature fields. 


Fig. 27 
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fig. 29 


Initial 1000 mb height (top) and temperature (bottom) fields for 
experiment GS valid at 12 GMT, 15 January 1979. 



Fig. 30 


Initial 1000 mb height (top) and temperature (bottom) fields for 
experiment GN /slid at 12 GMT, 15 January 1979. 




Fig. 31 



Initial 1000 mb height (top) and temperature (bottom) fields for 
experiment CS valid at 12 GMT, 15 January 1979. 
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Fig. 32 


initial 1000 mb height (top) and temperature (bottom) fields for 
experiment CN valid at 12 GMT, 15 January 1979. 
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The forecast experiments were all conducted using the 0.5° resolution 
version of the CMM described In Section 4. Forecasts were run out to 24 hours 
beginning at 12 GMT, 15 January 1979. Boundary tendencies for the 0.5° model 
were provided every two hours by the 1.0° resolution version of the CMM as 
discussed in Section 4. The resulting 12-hour forecasts for 1000 ab height 
and temperature for experiments VS, VN, GS, GN, CS, and CN are shown in Figs. 
34-39. The only real difference that is discernible from examination of these 
figures is the slightly stronger temperature gradient off the southeast coast 
for the assimilation cases. Subjective comparison of these fields is quite 
di f f i cul t . 

In order to objectively compare the 12- and 24-hour forecast fields 
produced in these experiments we shall compute root mean square errors for the 
differences between upper-air observations of height at 850 mb and 500 mb and 
the heights of the forecast fields interpolated to the locations of the obser- 
vations. We shall also use the FGGE Level Ill-b fields valid at the forecast 
times to compute SI scores (Teweles and Wobus, 1954) for the 850 mb and 500 mb 
heights. Since there are only about thirty upper air observations available 
over our domain spaced more widely than the surface observations shown in Fig. 
3, the root mean square errors we obtain are hardly conclusive. Similarly, we 
saw in Fig. 33 how smooth the FGGE Level Ill-b verifications fields can be. 
This is a limitation on the use of the SI scores. However, these data provide 
the only truly objective means to compare the forecasts. Table 2 summarizes 
the results of the statistical computations. At 850 mb for the 12-hour 
forecasts we can see that the. only noticeable result is that both VS and VN 
possess higher skill levels (lower SI scores) than the Cressman techniques. 



Fig. 34 


Twelve-hour forecast 1000 mb height (top) and temperature (bottom) 
fields for experiment VS valid at 00 GMT, 16 January 1979. 



Fig. 35 





Fig, 37 Twelve-hour forecast 1000 mb height (top) and temperature (bottom) 
fields for experiment GN valid at 00 GMT, 16 January 1979. 






Pig. 38 


Twelve-hour forecast. 1000 mb height (top) and temperature (bottom) 
fields for experiment CS valid at 00 GMT, 16 January 1979. 
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Fig. 39 


Twelve-hour forecast 1000 mb height (top) and temperature (bottom) 
fields for experiment CN valid at 00 GMT, 16 January 1979. 
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85£Hrib Height 

Forecast 


12— hour 

24- 

hour 


RMSE 

SI 

RMSE 

SI 

VS 

11,5 

26.9 

10.8 

34.8 

VN 

11.6 

26.9 

10.9 

35.2 

CS 

11,7 

30.8 

12.3 

38.6 

CN 

11.7 

31.0 

12.3 

38.8 

GS 

12.3 

32.4 

11.9 

38.1 

GN 

12.5 

32.2 

11.5 

38.4 



500-mb Height 

Forecast 



12- 

hour 

24-hour 


RMSE 

SI 

RMSE 

SI 

VS 

18.4 

13.3 

22.4 

17.3 

VN 

18.5 

13.5 

22.4 

17.2 

CS 

17,1 

13.3 

24.5 

18.1 

CN 

17.0 

13.3 

24.4 

18.0 

GS 

17,5 

33.2 

23.6 

17.8 

GN 

17,4 

12.9 

23.2 

17.6 


Table 2. 


Root mean square errors and SI scores for the forecast experiments. 



Difference fields between lowest 12-hour forecast model level 
temperature and temperatures interpolated from the FGGE Level Ill-b 
data set valid at 00 GMT, 16 January 1979, for experiments VS (top) 
and VN (bottom). 
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Fig. 40 
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After 24 hourn both the RMSE's and SI scores for VS and VN are lower than for 
the other experiments. At 500 mb there does not appear to be any significant 
differences among the SI scores or KMSE's for the different experiments 
although for 24 hours VS and VN appear to have slightly lower values. In no 
case do we see any indication of Improved skill or lowered RMSE's due to 
assimilation of satellite temperature data. However, due to the limitations 
of these objective comparisons, this Is not surprising. We saw before that the 
most noticeable effect of aosimilation was the enhancement of lower level 
temperature gradients off the southeast coast. From looking at Fig. 33 we 
know that this type of temperature gradient will not be reflected in the FGGE 
Level Ill-b data and hence the SI scores. Since the upper air stations are 
predominantly located on the land we cannot expect to see any reflection of 
this in the RMSE’s also. 

In order to more closely see the effects of satellite data assimilation 
we compared the differences between the temperatures for the lowest model 
level and for the FGGE Level Ill-b data interpolated to that level. In Fig. 
40 we can see these difference fields for the 12-hour VS and VN forecasts 
valid at 00 GMT, 16 January 1979. Similar resultB were seen for the other 
experiments. This figure points out that the only striking difference between 
these forecasts is the enhanced temperature gradient shown by VS off the 
southeast coast. This is a significant difference whose impact would not be 
detected by either of the objective statistical comparison techniques. 
However, given tiut the variational temperature analysis scheme produced 
enhanced low level temperature gradients along the southeast coast, we believe 
that this figure Illustrates that this information has been effectively 
assimilated into the CMM. 
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6. CONCLUSIONS AND RECOMMENDATIONS 

The satellite data available for this study were brightness temperatures 
from one orbit for the TIROS-N HIRS and MSU channels at 0850 GMT, 15 January 
1979. Upon examination of the data we concluded that It contained useful 
temperature Information In cloud-free areas for the very lowest levels of the 
atmosphere. The data were biased but were found to possess high resolution 
temperature gradient Information. Examination of the FGGE Level Il-b and 
Ill-b data we had available for that time showed that this was Information 
that could not be obtained from conventional sources. 

In order to combine the horizontal temperature gradient Information 
derived at high resolution from satellite with the absolute temperature 
Information obtained from conventional sources, we developed a variational 
temperature analysis scheme. This scheme was tested upon the data available 
and was found to be success. h'. at combining the two types of information. 
Using the scheme we were able to produce temperature fields which possessed 
the gradient information contained In the satellite observations and at the 
same time were unbiased and possessed relatively small RMSE’s with respect to 
the conventional observations of temperature. 

To test the Impact of the satellite data upon mesoscale model forecasts 
we conducted several forecast experiments utilizing the CIMMS Mesoscale 
Model. Since the domain of a mesoscale model is typically covered by one pass 
of a polar-orbiting satellite it wa6 natural to treat the data assimilation 
problem as an initialization problem. Three different model initialization 
procedures were tested with and without the use of satellite data. The first 
was a three-dimensional., multivariate variational objective analysis scheme; 


- 69 - 


the second was a two-dimensional, univariate Cressman scheme; and the third 
ws 8 identical to the second except that the resulting mass and momentum 
variables were coupled using a variational static balancing scheme. 

First, we examined the initialization procedures' effectiveness at 
providing balanced fields for the forecast model atart-up. Root mean squares 
of pressure tendency over the model domain were computed each time step as a 
measure of the initial imbalance and resulting model adjustment for forecasts 
initialized with each scheme. It was found that the variational scheme and 
the statically balanced Crewman scheme produced fields which permitted the 
model forecasts to begin with much less noise than the univariate Cressman 
scheme. Over the course of the model forecasts It was found that the noise 
level was less overall when the multivariate, variational procedure was used. 

Twenty-four hour forecasts were run for each of the three ini tiallzstion 
procedures with and without the use of the low level satellite temperature 
data. Objective statistical scores (RMSE's from conventional observations and 
SI scores from FGGE Level Ill-b fields) were computed and it was found that 
there was no significant impact upon any of the forecasts due to satellite 
data assimilation. However, it was found that the forecasts produced using 
the multivariate, variational initialize! ton procedure were significantly 
better than those produced using the two Cressman schemes. Subjective 
examination of the forecast fields revealed , however, that the satellite data 
assimilation did have an impact but that it could not be measured using 
conventional observations or analyses derived from conventional observations. 
All of the assimilation forecasts displayed the advection of the enhanced low- 
level temperature gradient from the southeast U.S. coast out into the 
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Atlantlc. Unfortunately, this feature, which vas the major difference between 
the assimilation and non-assimilation forecasts, occurre. an area with 

virtually no observational coverage. These forecasts did Bhow that the 

satellite-derived temperature information was indeed successfully assimilated 
into the forecast model. 

Based on the results of this study we would make the following 
recommendations for future research in this area. First, the brightness 
temperatures from the individual satellite channels were found to b^ useful 
only at the lowest levels of the atmosphere and for cloud-free conditions, A 
study similar to this one should be conducted using the highest horizontal 
resolution satellite-derived sounding data available from the TIROS-N 
satellite. The variational temperature analysis procedure was shown to be 
successful at combining the information obtained from satellite observations 
with that obtained from conventional observations. The impact of this was 
lessened, however, by the limitations of the satellite data available. In 
order to fully test its utility experiments must be conducted using a better 
collection of satellite observations. Finally, satellite-derived moisture 
profiles have been found to be similar to temperature profiles in that 
horizontal gradient patterns have been found to be reasonable and consistent 
with the meteorological processes involved (Smith, et al., 1981; Hayden et 
al«, 1981). An analysis scheme similar to the variational temperature 
analysis scheme should be developed and applied to satellite-derived moisture 
profiles and conventional observations. Forecast experiments similar to those 
described in this report should then be conducted to assess the impact of 


satellite moisture assimilation. 
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